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ABSTRACT 

We investigated wave amplification through the Electron-Cyclotron-Maser mechanism. 

We calculated absorption and emission coefficients without any approximations, 
also taking into account absorption by the ambient thermal plasma. A power-law en- 
ergy distribution for the fast electrons was used, as is indicated by X-ray and microwave 
observations. 

We developed a model for the saturation length and amplification ratio of the 
maser, scanned a large parameter space, and calculated the absorption and emission 
coefficients for every frequency and angle. 

Previous studies concluded that the unobservable Z-mode dominates in the u p w 
ub region, and that millisecond spikes are produced in the region Vp/vB < 0.25. We find 
that the observable O-mode and X-mode can produce emission in the 0.8 < Vp/vB < 2 
region, which is expected at the footpoints of a flaring magnetic loop. 

The important criterion for observability is the saturation length and not the 
growth rate, as was assumed previously, and even when the Z-mode is the most strongly 
amplified, less strongly amplified O-mode or X-mode are still intense enough to be 
observed. 

The brightness temperature computed with our model for the saturation length 
is found to be of order 10 16 K and higher. The emission is usually at a frequency 
of 2.06^3, and at angles 30-60 degrees to the magnetic field. The rise time of the 
amplified emission to maximum is a few tenths of a millisecond to a few milliseconds, 
and the emission persists for as long as new fast electrons arrive into the maser region. 

Key words: Sun : flares - Sun : radio radiation - masers - radiation mechanisms : 
non-thermal - instabilities - methods : numerical 
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1 INTRODUCTION 

Electromagnetic emission generated by relativistic electrons 
moving in a magnetic field is referred to as synchrotron 
radiation if the electrons are ultra-relativistic, and as gy- 
rosynchrotron if the electrons are mildly relativistic. In the 
present paper we deal only with gyrosynchro tron r adiation, 
using the formalism developed by Ramaty ( 1969), except 



that we also include the longitudinal component of the ra- 
diation which is important at low multiples of the cyclotron 
frequency. We also take in to acco unt the correction pointed 
out by Trulsen and Fejer ( 1970 ). This formalism has been 
tested by other workers (e.g. Holman & Benka 1992) and 
was found to be accurate. Gyrosynchrotron absorption can 
become negative and this amplification process is referred 
to as the Electron-Cyclotron-Maser (ECM). 

The observation of millisecond microwave spikes from 
the Sun (reviewed by Benz 1986) has been interpreted as 
gyrosynchrotron maser emission since the late 1970's (Wu 
& Lee 1979; Holman, Eichler & Kundu 1980), and it is now 
the most widely accepted explanation for the spike emission 



flBcnz 199$ . Spikes have been observed in the frequency 
range 0.3-8 GHz, exhibit a duration of a few milliseconds to 
a hundred milliseconds, a narrow bandwidth (0.5-2 per cent 
or a fe w MHz), and un til recently have not been resolved in 
angle (Benz et.al. 1996). By assuming that the homogeneity 
scale of the magnetic field is of order 10 7 cm, and that this is 
the size of the source, it is possible to deduce a characteristic 
brightness temperature of 10 15-18 K . 

In order that ECM appears the electron number distri- 
bution must either increase towards higher energies or an 
anisotropy in pitch angle must exist, as in a loss cone dis- 
tribution. A loss cone forms naturally where the magnetic 
field increases, as electrons with small pitch angles continue 
and electrons with large pitch angles mirror. At the foot- 
points of flare loops the density becomes high enough that 
electrons continuing into the chromosphere lose their en- 
ergy and do not reflect back upwards, and thus a loss-cone 
distribution appears (Aschwanden & Benz 1988; Kuncic & 
Robinson 1993). We therefore expect maser emission from 
the region of the footpoints. 

The study of the ECM mechanism is of interest because 
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of the radio emission, which serves as a probe into condi- 
tions i n the flare, and also because of the energetic particle 



(Aschwanden & Benz 1985), while in others the fun damental 



Ordinary Mo de (OM) is dominant for v p )vb < 0.5 (Winglee 



precipitation which produces X-rays, and for local heating 
of coronal material. The ECM reduces the trapping time 
and increases the flux of precipitating particles by causing a 
reduct ion in the p itch angle of electrons within a short dis- 
tance (Benz 199J ) . Electrons which otherwise would be mir- 



rored back and forth are thus precipitated into the chromo- 
sphere after one reflection. The maser thus transfers energy 
from the fast particle distribution to the ambient plama, and 
serves as a heating mechanism. 

Several papers have explored the possibilities of ECM. 
Melrose and Dulk ( 1982) have used a semi-relativistic ap- 



proximation to derive approximate formulae for the fre- 
quency of the la rgest growth rate, and for the growth rate. 
Winglee (| 198S|) approximated the effects of the ambient 
plasma temperature on the maser emission. Aschwanden 
(Aschwanden 1990a,b) followed the diffusion of electrons 
into the loss-cone as a result of the maser emission, and 
therefore the self-consistent closure of the loss-cone as ki- 
netic energy i s transf ormed to radiation through the maser. 
Aschwanden ( 1990b) concluded that the millisecond spikes 
are produc ed in the region of v p /vb < 0.25. Kuncic and 
Robinson ( 1993) performed ray-tracing in a model loop 
with dipole field, and concluded that maser emiss ion can 
escape from lower levels. Fleishman and Yastrebov ( 1994a) 
present results which show simultaneous amplification in two 
or more frequencies. 

An important difference in our work is that all previ- 
ous studies concentrated on the growth rate of the amplified 
wave, and defined the dominant mode as the particular wave 
- frequency v and direction to the magnetic field 9 - with 
the largest growth rate. However, the important criterion 
for long-term emission is the saturation length, which is the 
distance beyond which the m aser ceases to amplify, and not 
the growth rate (Benz 1993). We compute the absorption 



coefficient and, using our model, the saturation length, and 
determine which wave will be observed, and at what inten- 
sity. 

Another difference is that most previous studies used 
soft energy distributions for the fast electrons, such a s the 
Dory, Guest, Harris distribution (Winglee & Dulk 1986) or a 



hot thermal component distribution (Aschwanden 1990a,b). 
However, observations of Microwaves (MW) and X-rays 
from solar flares indicate a power law energy spectrum with 
energy up to 10 MeV and even higher, and with a power 
index typically S — 3, and in the range S — 2 — 6 (Ramaty 
1969; Zirin 1989; Fleishman & Yastrebov 1994b for other 
references). We therefore chose to use a power law distri- 
bution function for our fast non-thermal electrons. We also 



use an ideali zed one sided loss-cone distribution ( Melrose 
& Dulk 1982) while other studies use a one sided sin (<f>) 
distribution (Aschwanden 1990a, b) or a two-sided gaussian 
(Fleishman & Yastrebov 1994a,b) pitch angle distribution. 

The ratio v p /ub, where v p = y n e /m/-n e is the plasma 
frequency and vb = eB / (rwyc) / (2n) the cyclotron fre- 
quency, is an important parameter for ECM emission. It 
is used to distinguish regions where different magneto-ionic 
modes dominate. From the previous studies we conclude 
that the dominant mode depends on the distribution func- 
tion of the fast electrons. In some studies the fundamen- 
tal eXtraOrdinary (XO) mode is dominant for v p /vb < 0.3 



Dulk 1986). In some the fundamental OM is dominant for 
0.3 < v p )vb < 1 (Aschwanden & Benz 1986), and in others 



t he second harmonic XO is dominant for 0.5 < v v /ub < 1.5 
( [Winglee fc Dulk 1986~| ). The unobservable Z-Mode ( ZM) 
is dominant for v p /ub > 1.5 (Winglee fc Dulk 1986 ), or 
even for v p /vb > 0.5 (Fleishman & Yastrebov 19941). In 



our work the ZM is the most highly amplified mode (i.e. 
with the largest in absolute value negative absorption) for 
the range 0.5 < v p /i'b < 1.1 and for 1.8 < v p /vb, similar 
t o what was found by others w ho have used a power law 
( Fleishman fc Yastrebov 1994b ). 

The high amplification of the ZM would seem to in- 
dicate that this unobservable mode is the dominant mode, 
and would therefore preclude the maser mechanism as the 
producer of MW spikes, except in the narrow range 1.1 < 
v p Ivb < 1.8. We show that for the range v p /vb < 2 it is 
possible for the OM or the XO mode to be amplified enough 
that they become observable. Our absorption-coefficient ap- 
proach allows a simple procedure to estimate the satura- 
tion length, which is difficult to derive from growth rate 
calculations, and is the important parameter in deciding 
which mode dominates (Benz 1993). Our derived satura- 
tion lengths are consistent with the assumed homogeneity of 
the magnetic field and density. We therefore conclude that 
the Electron-Cyclotron-Maser mechanism can produce mi- 
crowave spikes. 

In section 2 we give the theory and equations for gy- 
rosynchrotron emission and wave propagation through cold- 
plasma. In section 3 we show that the gyrosynchrotron is the 
important process in absorption as well as emission, for the 
frequencies and plasma parameters we assume. Section 4 de- 
scribes our estimate of the energy available to the loss-cone, 
and a derivation of the saturation length is given. Section 
5 describes our calculations and the parameters we used. In 
section 6 we discuss the observability of the maser emission 
under different conditions, in section 7 we give a compar- 
ison with previous studies, and in section 8 we summarize 
our results. 



2 GYROSYNCHROTRON ABSORPTION 

Solar flares are known to be a complex phenomena. The so- 
lar flare plasma is composed of both an ambient thermal 
plasma, and a population of fast particles. It is customary 
to simplify and assume that any waves which are excited 
are a phenomena of the ambient plasma alone. The proper- 
ties of the ambient plasma therefore determine which waves 
are excited and how the waves propagate. However, the fast 
particles can absorb waves present in the ambient plasma, 
and also emit in these waves' frequencies. We used the cold 
plasma approximation to describe the ambient plasma, and 
assumed the w aves which result from magnetoionic theory 
( Melrose 1989 ). For a given electron distribution f(p) and 
wave (v, 8) it is then possible to compute the emission and 
absorption coefficients, which describe the interaction of the 
wave with the distribution. 

The cold plasma parameters which determine the prop- 
agation are the refraction index and the components of the 
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polarization vector (Melrose 1989) 



L a (k)k + T a {k)t + ic 



(1) 



where k is the wave vector, k a unit vector in the k 
direciton, c a unit vector perpendicular to k and to the 
magnetic field B, and t is a unit vector perpendicular to 
k, in the plane of k and B. Also T CT and L a are functions 
depending on k, and a is '+' (OM) or '-' (XO mode). The 
transverse component of the polarization is given by T a and 
the longitudinal component is described by L a ■ We assumed 
B to be in the z direction, and k to be in the (x, z) plane 
with angle 9 to the magnetic field. With these definitions 
c = y , k — (sin(#), 0, cos(0)), and t — (cos(#), 0, — sin(0)). 

The cold plasma modes are the OM, which exists for 



v > v„ 



the XO which exists for 

v > v x = vb/2 + \Jvl + v 2 /4 



(2) 



(3) 



and the low frequency branch of the XO, called here the ZM, 
which exists for 



vb 



< v < V0.5 x 



+ "D + [("p + "D* - CO S 2(0)]O.5 

and is bound from above by the upper-hybrid frequency 

(4) 



Vub = \/ U B + U P 

The Whistler mode exists for 



v < 



(5) 



and is a low frequency mode which we ignore (Melrose 1989) 
The equations for the cold-plasma parameters are (Ra- 
maty 1969; Melrose 1989) 



[2z/(z/p - u) cos(6>)] x 
[ — v 2 vb sin 2 (9) ± 



4,,2 



v 2 ub sm{9) 



\9)+hj 2 {v 2 ~v 2 ) 2 
vT ± 



UUp vT± — ub cos(S) 



(6) 



(7) 



And the refraction index n rc f, also referred to as n±, is 
given by 

n\ = l + [2vl{vl-v 2 )} X (8) 
[ — 2v 2 (v^ — v 2 ) — v 2 v\ sin 2 (S) ± 

sin 4 (0) + Av*vl{vl - v 2 ) 2 cos 2 (6») ] _1 

With the symbols being v p the plasma frequency, vb 
the cyclotron frequency, and v the emitted frequency. 

The cold plasma approxi mation is valid only if the two 
following criteria are are true ( fMelrose 1989] ) 



n IC {V cos(9) 
and 
vb 

n ro { v sin(#) 



>> 



k B T 

m e c 2 



0.001683 



T 



10 7 K 



>> 



k B T 



0.001683 



10 7 K 



(9) 



(10) 



Where we have used k\\ = n re fW cos(9)/c and wthormai = 
yj feT/m e to get the above from Melrose's equations. The 
s multiplier in the first criterion is an integer such that sub 
is the closest cyclotron harmonic to v. 

Our rule for when the cold plasma approximation ceases 
to be valid is when the left side of the criteria is less then ten 
times larger than the right side. The cold plasma approxi- 
mation fails near the resonances of the cold plasma modes at 
the cyclotron harmonics, and for the low-frequency Z-mode 
when the frequency of the waves is near the upper-hybrid 
frequency. 

Most authors calc ulate th e growth rate, but our deriva- 
tion follows Ramaty's ( 1969) in calculating the absorption 
coefficient which is (Ramaty 1969; Aschwanden 1990a): 



B \cos(9)\ v n £ ± 



' j u (7) 
d7 x 



(11) 



(3 sin ( 



I J3 (x s ) -I J s (x a 

n± 



+ 



I cot [ff) P COS (<j) B ) 



-Pj 2 d u (7) 



sin (9) 



u (7) d7 /?7 2 
n±j3 cos (9) — cos ( 



Js (X s ) 



dg(4>) 



7/3 2 sin (4> s 



Where we corrected Ramaty's ( 1969 ) formula by delet- 
ing a factor 1 + d\n{n re f) /d\n{v), and changed the minus 
sign multiplying the J' s term because our T± have the oppo- 
site sign to Ramaty's ag±. We also added the longitudinal 
contribution which is important for frequencies close to the 
cyclotron frequency. J s are Bessel functions, their deriva- 
tives, with the argument x s = (^/i^B)n±7/3sin(#) sin(0 s ). 
The limits of the summation are si, 2 = ^ [l±n fe f/3 cos(#)], 
and ± designates the mode ('+' for OM and '-' for XO). The 
angle <j> s comes from the resonance condition (Ramaty 1969; 
Melrose 1989) 



— — = f[l — n rc f/3 cos(#) cos((j>)] 



and is derived by 
1 



(12) 



(13) 



n± cos(9)P 

The distribution function of the electrons was assumed 
to be separable into the form 



f(p)d 3 p = 27TM(7)p(^)d7d(cos((/))) 



(14) 



Where <j>, the pitch angle, is the angle of the velocity to the 
magnetic field, and 7 = (1 — /9 2 ) -0 ' 5 the Lorentz factor. 

It is important to note that the equation for the absorp- 
tion coefficient has the same properties as the equation for 
the growth rate since the growth rate is 



r± 



-fe± ■ v g 



(15) 



where v s is the group velocity given by (Melrose 1989) 
du) c 



dk g(""rcf) 



(16) 



1 + T 2 

As in all previous studies, we neglected the angle between 
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the group velocity and the wave vector. The group velocity 
in our study is therefore assumed to be in direction k, and 
of magnitude 

(17) 



In this work we assumed the power-law electrons have 
a power-law dependence on kinetic energy 



w(t) = rihot ■ {5 - 1) • (7mm - 1) S 1 ■ (7 - 1) S 



(18) 



where rihot is the number density of the fast electrons, and 
the power law starts from a low-energy cut-off 7 m i n and ends 
at a high-energy cut-off 7 max . The distribution is normalized 
to n ho t when we neglect (7 m i n — l)/(7max — 1). For the am- 



bient plasma we used a thermal distribution (Lang 1980) 



i(7)d 7 = n cold -^= ( ) y/y- 1 



3/2 



x e kT 



(7-1) 



d7(19) 



where n co id is the number density of the cold, ambient 
plasma, electrons. 

For the pitch angle de pendence in the fast d istribution 
we used an ideal loss-cone ( Vlelrose fc Dulk 1982 ) 



Const.. 

Const * cos ( a - Aa ~)- cos (<t>) 

g(<fi) = < ' ms(o-Ao)-cos(a)' 



0, 



cos(4>) < cos(a) 
cos(a) < cos(<^>) < 
cos(a — Aa) 
cos(a — Aa) < cos(0) 



This describes an isotropic distribution for pitch angles 
above some a, no electrons below another angle a — Aa, and 
a linear in cos((f>) decrease within the Aa interval. The nor- 
malizing constant in g(<j>) depends on the loss cone param- 
eters. For an idealized loss cone as above it is Constant = 
l/7r/(2 + cos(a) + cos(a — Aa)). With this definition the 
integral is normalized so that 2ir J g(if>)d(cos(f>) = 1. 

For the thermal plasma the distribution is isotropic and 
therefore g(<jj) = 1/4/n. 



3 OTHER PROCESSES 

There are other processes in a thermal plasma which have 
to be taken into acco unt. The fr ee- free absorption coefficient 
for T > 3 • 10 5 K is flLang 198C|) 



fc « = °- 01 ^( 24 - 57 + ln ?) 



(21) 



where n e is the density of the thermal electrons in cm , v 
the frequency of the absorbed radiation in Hz, and ks the 
absorption coefficient in cm" 1 . 

Our assumed number density is of order n e < 10 , 
our standard temperatures are of order T « 10 r , and the 
emission frequency for our standard values is of order v ~ 
1GHz. The free-free absorption for these values is fcg « 5 ■ 
10 -8 . For a standard number of fast electrons of order nhot ~ 
10 6 cm -3 the free- free absorption normalized to one electron 
is fcff/nhot < 5 ■ 10~ 14 , which is negligible compared to the 
gyrosynchrotron absorption of the dominant modes. 

The deflection time for a particle in a thermal plasma, 
with a particle velocity much hig her than th e thermal ve- 
locity but still non-relativistic, is (Benz 1995) 



t d = 3 -10- 20 ^ 

n e 



(22) 



(20) 



where v is the particle velocity, and ta is the time for a 
deflection of 90 degrees. 

Taking an electron with 7 = 1.02, which is the low 
energy cutoff of our fast electron distribution, the deflection 
time is td ~ 62 msec. The time for a deflection of only a few 
degrees is much shorter, but the timescale of saturation is a 
few milliseconds, and therefore the deflect ion should not be 
an important factor ( Aschwanden 1990b] ). Faster particles 
have even longer deflection times, and are even less affected. 



4 CONDITIONS FOR MASER AND 
SATURATION LENGTH 

Inspection of equation [ll] shows that the absorption for a 
distribution which is isotropic in pitch angle is negative only 
if the slope of 14(7) /j 2 //3 is positive. This is the equivalent of 
an increase of the fast particles with energy, a distribution 
which is not likely to occur. The other possibility for am- 
plification is when the distribution is not isotropic in pitch 
angle. In that case if the derivative of g(4>) is large enough, 
and there are solutions such that the pitch angle dependent 
term in equation [ll] is negative, then the absorption will be 
negative. One such distribution is the loss-cone distribution, 
which we use. The loss-cone distribution arises naturally in a 
solar flare where the magnetic field and the density increase 
with decreasing height from the corona to the chromosphere, 
and we expect a loss cone to appear near the foot-points of 
a flare. 

A numerical computation of equation [ll with a loss- 
cone (equation ^) and a power-law (equation ^) results in 
negative absorption coefficients for several frequencies and 
angles to the magnetic field. Waves of these frequencies, and 
propagating with these angles, are amplified, and this is the 
Electron- Cyclotron-Maser mechanism . 

The major factor determining which mode is the most 
important, or dominant, is k x L, or F x L/v s , where L is the 
length of amplification (Benz 1993). If L is determined by 
homogeneity length-scales in the loop then the mode with 
largest absolute value negative absorption is dominant. How- 
ever, if the maser stops operating after a distance smaller 
than the homogeneity scales of the loop, the amplification 
length is intrinsic to the ECM itself. Our goal was to derive 
the absorption coefficients, and the length of amplification 
L. 

The intrinsic length of the maser was estimated with 
a model for the evolution of the fast electron distribution. 
We assumed the electrons are accelerated for a short time 
at the loop top, so that a pulse of fast electrons with a du- 
ration of a few milliseconds to a few hundred milliseconds is 
produced. The fast electrons reach an abrupt increase in the 
magnetic field, and the electrons with large pitch angles are 
reflected back. We assumed that before reaching the mirror 
the fast electron distribution is isotropic in pitch angle, and 
that all the electrons below a critical pitch angle pass the 
field increase. These assumptions result in the ideal loss cone 
distribution described by equation |io| with a loss-cone open- 
ing angle of a. Our evolution model therefore began with an 
ideal loss-cone distribution over some volume, and at time 
t = the maser is "turned on". We then assumed that the 
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anisotropic loss-cone distribution relaxes into a new distri- 
bution via the maser, and for simplicity we assumed the new 
distribution is also an ideal loss cone. The new distribution 
has a smaller opening angle at, and a power law in energy 
of the same power index as the original distribution. If the 
dominant process which causes the relaxation is the maser, 
the difference in energy between the initial loss-cone distri- 
bution and the final distribution is the energy available to 
the maser. 

A useful measure of the amplification process is the sat- 
uration length, i sat . We define the saturation length as the 
distance a wave can be amplified before the loss-cone dis- 
tribution is relaxed so that the loss cone disappears, or the 
particular wave is no longer amplified. Complete relaxation 
for a volume element occurs when the energy in the ele- 
ment is reduced from the anisotropic distribution level to the 
isotropic distribution level. Our method for calculating the 
saturation length assumes that the greatest possible ampli- 
fication is when the entire energy difference between initial 
and final loss-cone configurations is turned into maser radi- 
ation. We estimate this energy difference by assuming that 
the relaxation process leaves the number of fast electrons 
and the total parallel momentum unchanged. The length of 
amplification can not be larger than the smallest saturation 
length, and therefore if i sa t is smaller than the loop length- 
scales it determines the dominant mode. The mode with the 
smallest saturation length is therefore the dominant mode. 

We followed the gradual closing of the loss cone in detail 
for some test cases, and our results justify an all-or-nothing 
approximation. In this approximation we assume that the 
distribution remains unchanged until the energy radiated 
by the maser is equivalent to some per cent of the kinetic 
energy, and then the relaxation to the lower energy configu- 
ration is instantaneous. When the closing of the loss-cone is 
followed it is found that the loss-cone opening angle is usu- 
ally changed only by a few degrees before the absorption co- 
efficient for the mode which had the largest amplification be- 
comes positive. However, this happens after a distance which 
is only larger by 20 to a 100 per cent of the Z sa t calculated 
with the all-or-nothing approximation, when we assume the 
energy turned into the maser is 10 per cent of the kinetic en- 
ergy. The energy difference with the gradual closing method 
is usually 7-10 per cent of the total kinetic energy. Our all- 
or-nothing method is simple and gives a transparent relation 
between the absorption and the saturation length. 

Our calculations are predicated upon the assumptions 
that the magnetic field strength B, the direction to the ob- 
server 6, and the densities ra co id, rihot are constant. Figures 
^, ^| and [B| show contours of negative absorption coefficients 
around the area of maximum negative absorption. From 
these figures it is evident that a change in direction of the 
magnetic field by even 1° removes a propagating wave from 
the peak of negative absorption. Likewise a change of 1 per 
cent in the magnetic field makes the normalized frequency 
v/vb of a wave change by about 1 per cent, and also removes 
the wave from the peak of amplification. Estimates of the 
scale lengths for the magnetic field are usually such that 
chan ges of 1 per cent occur over distances of about L — 10 7 
cm (Aschwanden & Benz 198S ) . Therefore, for saturation 
lengths much larger than 10 cm, we expect the amplifica- 
tion ceases before the particular wave reaches its saturation 



length, and the loss-cone either persists through the loop, 
or perhaps relaxes by some other means. 



4.1 Energy in the Loss-Cone 

The total momentum parallel to the magnetic field for a one- 
sided loss cone distribution with constant density in angle 
up to cos(<f>) = cos(a) and no electrons for higher cos(<j!>), 
and with power law in kinetic energy is 



(23) 



= 71-/2 



d7d(cos(<^))m7/3ccos(0)27ru(7)g(<^) 



mc cos 2 (a) 



2tt(1 + cos(a)) 2 
Where the dependence on energy is included in the integral 



W 7 = n ho t(<5-l)(7» in -l)' 5 - 1 / ^ 7 2-l( 7 -l)-*d 7 (24) 

and 7 m in is the lowest energy, 7 max the highest energy, of 
the fast electrons. 

The dependence on the opening angle of the loss-cone 
a is included in the first part of the formula for Pn, where 
the integral over cj> was done on the positive side of the dis- 
tribution, 4> < 90, and where we assume the transition zone 
from isotropic distribution to the empty region is negligible. 

The energy dependent integral W n is especially simple 
for the case of a power law in energy with 5 = 3 and is 



2rlhot ,., \3/2/ i\l/2 

VV 7 = — - — (1 + 7min) (7mm - 1) 



(25) 



where we have neglected, as is usually done, terms of order 
(7min — l)/(7max — 1). We use 7 m in = 1.02 and 7 max > 3 in 
our calculations. 

We also need the equation for the kinetic energy density 
in a power law distribution which is 



= 7 — -mc rihot (7n 



1) erg cm 



(26) 



where rihot is the number density of fast electrons. 

We demand conservation of parallel momentum when 
the distribution changes from a loss-cone with open- 
ing angle a and 7 m m to a distribution with j m i n r and 
at, so that P|| (cos(a), 7 m i n ) = P\\ (cos(a/), 7 m in/). From 
the parallel momentum conservation condition, and given 



cos(a), 7 m in, cos(q/), we calculate 



In 



From 7min' and 



equation B6j we calculate the amount of energy lost in the 
transition from a loss cone with a to a loss cone with opening 
angle at. The transition to a lower cut off energy is physi- 
cally reasonable since we expect electrons with low 7 to emit 
in the low frequencies more than electrons with high 7, and 
when these electrons lose energy they move to smaller 7-es. 

As an example we took an initial distribution with 
cos(o:) = 0.8, 7min = 1.02, which is similar to our standard 
loss-cone. From equations |2^ to ^ the isotropic distribution 
(cos(q/) = 1) which has the same total parallel momentum 
is a distribution with 7min' = 

1.01026. Therefore, about fifty 
per cent of the energy was lost in the transition from a loss- 
cone with cos(q)=0.8 and (7 m i n — 1) = 0.02, to an isotropic 
distribution with cos(q/) = 1 and (-/min r — 1) = 0.01026. An 
estimate of fifty per cent is two orders of magnitude larger 
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than that found by Mackinnon, Vlahos &Vilmer ( 1992) who 



estimated the energy as 1 per cent of total kinetic energy, 
and one order of magni tude larger than that given by White, 
Melrose & Dulk ( 1986). Howeve r, this estimate is about the 
same as quoted by White et.al. ( 1986) as resulting from an- 
alytic arguments. 

We also expect that a particular wave is no longer am- 
plified once the loss cone characteristics have changed even 
slightly. Therefore we should have calculated the energy loss 
when the loss cone angle changed from cos(a) = 0.8 to, say, 
cos (a/) = 0.82. The energy loss for this case turns out to 
be only seven (7) per cent. Since our purpose was merely 
to get an estimate for the energy available to the maser, we 
can take as a reasonable estimate ten per cent of the kinetic 
energy. The exact number has very little effect on the sat- 
uration length, but the brightness temperature is directly 
proportional to the energy (see equation S). 



4.2 Saturation Length 

To calculate the saturation length from the energy estimated 
to be available to the maser we assume that the maser is 
dominated by one mode {v,0). We assume that the entire 
energy is concentrated in this one frequency and direction, 
with a typical angular width according to our numerical 
computations of AO. = 0.1 it, and a typical frequency width 
of Au = 0.015^ B . 

To decide which is the most important mode we used 
the typical frequency and angular widths, calculated the sat- 
uration length for every mode as if that mode was the domi- 
nant one, and then compared saturation lengths. Usually the 
mode with the largest absolute value absorption coefficient 
also has the smallest saturation length, and is therefore the 
dominant mode. 

A justification of this approach is seen for example in 
Fig. ^[ The figure shows contours of equal absorption co- 
efficient for the case v p /vb = 1,T — 5 • 10 6 , and the OM. 
The absorption coefficients are negative, and therefore the 
figure shows equal amplification contours. The maximum 
amplification in Fig. |^ is in the region 0.795 < cos(#) < 0.8 
and 1.04 < vjve, < 1.045. If the magnetic field strength 
decreases by one per cent the frequency of the wave moves 
from approximately 1.04z/b to approximately 1.05z/b, and 
the amplification is decreased by approximately half. Like- 
wise, if the inclination of the field changes by one degree, 
the amplification is also reduced by approximately half. It 
is therefore justified to use narrow frequency and angular 
widths. 

The simplest method of calculating the saturation 
length is the all-or-nothing method. The amplification of 
any mode is calculated from the length over which there 
is negative absorption, which is assumed to be arbitrary at 
first. On the other hand the intensity produced from a given 
volume must be limited by the total energy in the volume. 
The length for which the loss cone persists is calculated self- 
consistently by equating the energy in the wave, and the 
energy in the volume. 

The total power of the maser per unit area is 



P = I(v,cos(6))AQAv 



(27) 



The energy increases over a volume equal to the path of 



amplification times unit area, and therefore the total energy 
available to the maser is 



0.1 x E k x Z sat [erg cm ] 



(28) 



where we used section 



1.1 



and our estimate of 10 per cent of 
the total kinetic energy as the energy available to the maser, 
and Zsat is the length of amplification - which is unknown as 
yet. 

Assuming that the maser lasts for a given time At, and 
that the intensity of the emission is constant for that time, 
the equation which has to be solved is 

PAt = 0.1x£ k xl sat (29) 

By using the radiation transfer equation for negative 
absorption 



J_ \k\u 

Ifcl 



(30) 



where we have neglected 1 compared to the exponent, and 
assumed that j, k are constant everywhere, we have 



^AQAvAt - 



J r ^hot\ k \ l B 

0.1 x 2 x 8.2 • 10~ 7 x n hot x (7 min - 1) X l a 



(31) 



where we used equation |26j with 3 = 3, and also introduced 
normalized (to 1 particle) absorption and emission coeffi- 
cients. 

If we assume the maser lasts for a given time At the 
equation can be solved numerically for any At. However, if 
we assume At = i S at/w g) with v g given by equation then 
the solution is straightforward, and the result of the simple 
method is 



/sat = In 



0.1-E k -v g - | 
AflAu ■ j 



/(nhot |*|) 



(32) 



A more complicated method of deriving the saturation 
length is the differential approach. We assume a loss-cone 
distribution over some volume of space. We further assume 
that the maser is "turned on" at time t — 0. The intensity at 
time t arriving into a volume element is given by equation 
|3C| with x = v g t, since a wave starting at time t = has 
advanced this distance by time t. A volume element adds to 
the intensity a fractional intensity 



dl = (j - kl)dx 



(33) 



and again neglecting 1 compared to the exponent, the in- 
tensity increase for a volume element is 



|U|fc|/«je' fc ^ 



(34) 



The power in erg s _1 drawn from a unit volume at time t 
after the maser was turned on, assuming that k and j remain 
independent of t and the same everywhere, is therefore: 



AP(t) 



dfl / duje 



k \vgt 



(35) 



The energy drawn from from a unit volume from time t\ to 
time ti is 

AE = j AP(t)di« / d« / dv-i— (e |fc| ^ t2 -e |ft| " Etl )(36) 

! H J J 1*1% 
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Since energy is being drawn from the volume element 
the distribution function changes, and here we use section 



4.1. We assume that the distribution changes so that it re- 
mains a power-law with index 8 in energy, and a loss-cone. 
The new distribution function, with less energy, has a new 
and smaller 7mm' and a new and larger cos(a/). From our 
calculated energy loss for some interval ti — t\ we write 
E' k = i?k — AE, and derive 7'. From the condition that 
the parallel momentum is constant P/u = Pu, and using 
7', we derive a' . We can then calculate the new absorption 
and emission coefficient, the new energy loss, and so on, 
and repeat the process iteratively. If we take t\ = and 
t2 = t = iaat/wg, and substitute for the integrals the con- 
stants AO. and Av, the equation for the saturation length 
we derive from the differential approach is identical to the 
equation of our simple method (equation ^) . 

We followed the gradual closing of the loss-cone itera- 
tively for a few test cases. We took into account the change 
of the absorption and emission coefficients for waves arriv- 
ing in later times, which passed through regions where the 
distribution has already changed. The results of this method 
are that about 3-7 per cent of the kinetic energy is turned 
into maser radiation, mostly in the dominant mode, that this 
happens after a time < 2 • r, and that the intensity of the 
maser is very near to the intensity after amplification with 
the original maximum absorption over length l sat = TV g . 
Thus the natural growth time of the maser to maximum 
is < 2 r and the length over which maser amplification is 
possible is < 2v s t, and our simple method is a very good 
approximation. 

Since the results of following the closing of the loss-cone 
and just assuming a fixed energy available to the maser are 
so close, we used for the rest of our calculations a fixed 10 
per cent energy availability. The results given in this work 
were calculated with this assumption. 



5 CALCULATIONS 

Our calculations were carried out within the frame of the 



cold- plasma and the magneto-ionic approximations (Melrose 



198S). We therefore disregard any negative absorption for 



frequencies which do not obey either of the two cold-plasma 
criteria ^ and |l^, where "much larger" in these criteria is 
taken to be at least 10 times larger. 

We calculated the absorption coefficient with equation 
[y], and the emission coeffi cient with the equation appro- 
priate for it ( Ramaty 1969 ). We performed the calculation 
using a power-law distribution in equation [n] to get the ab- 
sorption due to the fast electrons, and also using a thermal 
distribution in the equation to get the absorption due to 
the ambient plasma, and the absorptions were then added 
to derive the total absorption. The total emissivity was de- 
rived in the same manner. The calculations were carried out 
without any approximations to the Bessel functions, and the 
full sum over the s values was performed. The integral over 
the pitch angle (j> was performed using the fully relativistic 
resonance condition n2l and the distribution functions f(p) 
were normalized to the number density of the fast electrons 
for the power law, and the ambient plasma for the thermal 
distribution. The parameters used for the calculation were 



• 7min = 1.02 ; 

• 7max — 3, 

• cos(q) = 0.81 ; 

• cos(a — Aa) — 0.83 ; 

• B = 360 gauss ; 

• w C oid/»T-hot = 10 4 , unless otherwise indicated ; 

• The ambient density is taken to be n co id = 1.24 ■ 
10 10 (v p /vb) 2 cm -3 , which gives v p = 1 GHz for v p /vb = 1. 

We used 7 max = 3 because several tests we performed 
showed that electrons with higher kinetic energy do not con- 
tribute to emission in the frequencies of the maser. 

For every given pair of (v p /vb,T) we calculated the 
absorption and emission coefficients on a grid of cos(#) = 
0.02 - 0.92, with a spacing of A cos(0) = 0.02. For ev- 
ery cos(#) we calculated the coefficients for the frequency 
range from v = vb to v — 3vb, or higher if required, and 
with a changing resolution. The initial scanning resolution 
was changed from Av/vb = 10~ 4 for the small cos(6 l ) to 
AvjvB = 0.01 for the larger cos(6>). When negative absorp- 
tion coefficients were found the resolution was increased ten 
fold, and the calculation for the range where the absorp- 
tion is negative was performed with the higher resolution. 
In this way we are confident that no negative absorption in 
the (cos(6 l ), v) range was missed. Sometimes calculation was 
performed with better resolution in A cos(6>) or Av/vb if the 
absorption coefficients changed too much with the coarser 
resolutions. We performed calculations for the entire range 
Vp/vs = 0.1 — 3, but present only examples. 

Negative absorption appeared for the OM and XO near 
the harmonics of the cyclotron frequency, and for the ZM 
also near the upper hybrid frequency. Therefore it was usu- 
ally not necessary to search the whole range vb — 2vb and 
the whole range 2vb — 3^b to find the negative coefficients. 

Calculation of the saturation length was performed as 
described in the section Saturation Length, with the assump- 
tion of constant width in angle and frequency. The widths 
enter the saturation length only logarithmically, so the re- 
sult is only weakly dependent on them. Typical numbers are 
An = O.Itt and Av = 0.015f B . 

The tables summarizing our results are arranged ac- 
cording to Vp/vB- Every row gives the characteristics of the 
most important wave for the given mode and temperature, 
and for the v p /vb of the table. The waves are either those 
with the largest (in absolute value) absorption coefficient, 
or those with largest growth rate, for all 9 and v. The first 
column is the temperature used to calculate the absorption 
coefficient for the ambient plasma, in units of 10 6 K. The 
temperature affects only the thermal absorption and not the 
absorption calculated from the power law electrons. The sec- 
ond column is the frequency in units of the cyclotron fre- 
quency vb- The third column is the emission angle to the 
magnetic field, or the angle of the wave vector k (assumed 
to be in the (x, z) plane) to the z-axis ( B was assumed 
to be in the z direction). The fourth column gives the ab- 
sorption coefficient itself, normalized to a number density 
of the fast particles of 1, and given for a magnetic field of 
360 gauss which gives vb = 1 GHz. The absorption coef- 
ficient is therefore nhot times the k which is given in the 
table. The fifth column is the saturation length in units of 
10 6 cm, computed with our formula 



that nho 



1.24- 10 b {vp/vbY 



I. and the assumption 
The sixth column is 
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the group velocity, the seventh is the growth rate, and the 
last column gives the mode for which these parameters were 
computed. 

The contour figures show only negative absorption con- 
tours, multiplied by a factor of — 10 12 to make the numbers 
easier to display. These are examples showing only the re- 
gion of (cos(#), v/vb) plane around the largest in absolute 
magnitude negative absorption coefficient for some mode. 



6 OBSERVABILITY OF THE MASER 

A previous study which used a power law in energy for 
the fast electrons found that the ZM dominates for all 



Vp/vs > 0.5 (Fleishman & Yastrebov f994b), and therefore 
concluded that the emission must come from regions with 
smaller i/ p . Another study also concluded that the maser 
must come from a region of low density p p /ub < 0.25, be- 
cause for higher densities the dominant modes reach satura- 
t ion much faster tha n the observed durations of MW spikes 
( |Aschwanden 1990b| ). 

We found that for v^jvB < 1.1, and u^/pb > 1.8 the 
ZM has the largest in absolute value negative absorption (see 
Fig. hi), but has the largest growth rate only for v v )vb < 0.6 
(FigTbh. If the pulse of fast electrons has a much longer 
duration than the time scale of maser saturation the im- 
portant factor is not the growth rate, but the saturation 
length. At a time after the maser saturates, the mode with 
the smallest saturation length has closed the loss-cone be- 
yond its own saturation length, and therefore the smallest 
saturation length becomes the amplification length for all 
modes. Modes with longer saturation lengths can reach only 
a fraction of the full amplification possible to them, being 
amplified by a factor exp(\k\-L), with L the saturation length 
of the mode with the smallest saturation length, and not 
cxp(|fc| • igat), with Z sat the saturation length of the mode 
which is L < Z sat . Therefore what is required is to derive the 
saturation length for all modes, and find the mode with the 
smallest saturation length. Since the saturation length is in- 
versely proportional to the absorption coefficient, the mode 
with smallest saturation length is usually the mode with the 
largest (in absolute value) negative absorption coefficient. 

These findings present a problem in that we expect 
a loss-cone distribution to appear near the foot-points of 
flare loops where the magnetic field and the density increase 
sharply. The cyclotron frequency-plasma frequency ratio for 
the footpoints is probably in the range p^/vb ~ 0.5 — 3, if 
we take as typical w C oid ~ 10 10 - 10 11 cm fo r the density 



of ambient plasma (Aschwanden & Benz 1997), and B ~ a 



few times 100 gauss as typical for the magnetic field. These 
values are in the range where our calculations, and previous 
work, predict the unobservable ZM is the dominant mode. 
The high amplification of the ZM would therefore seem to 
preclude the ECM mechanism as the producer of observable 
MW spikes, except in the range 1.1 < v p /vb < 1.8. 

We have, however, to consider the ambient thermal 
plasma at the foot points, which we assume to be at a tem- 
perature of order 10 7 K. The thermal effect was not taken 
into account in most studies, and those which did - assumed 
much lower temperatures for the ambient plasma than we 
do. We use the range of temperatures inferred from SXR ob- 
servations, since we assume the model of a hot dense plasma, 



heated by the fast electrons, flowing upwards to fill the loop, 
and emissi on where the fast electrons intersect the evapora- 
tion front (Aschwanden & Benz 1997). 

The temperature has two effects on our results. The 
first is to increase the positive absorption near the harmon- 
ics of the cyclotron frequency. If the frequency where the 
fast electrons create negative absorption is close to a har- 
monic, the thermal plasma absorption decreases the am- 
plification, or even turns the negative coefficient positive. 
For the XO and OM this may eliminate the amplification 
at a low harmonic, but these modes appear in higher har- 
monics and usually have negative absorption for the higher 
harmonics as well. The ZM is limited to frequencies less 
than the upper-hybrid frequency, and therefore it is possi- 
ble that the introduction of high temperatures will eliminate 
the ZM negative absorption while leaving higher frequency 
OM and XO negative absorptions. The second, and more 
important, consideration high temperatures introduce is the 
cold plasma conditions. For a temperature of 10 6 K these 
criteria are valid for frequencies very close to the harmon- 
ics, and for very large refraction indices. However, the high 
temperature we assumed forced us to reject certain negative 
absorption results as not valid within the framework of our 
assumptions. We assumed that if there is any amplification 
for these cases it is orders of magnitude smaller than the cold 
plasma result, and therefore negligible. This is supported by 
the known effect of thermal corrections, which is to decr ease 
the refraction indice and the absorption ( Melrose 1989| ). 

Even if the thermal positive contribution leaves a neg- 
ative ZM absorption which is valid for the cold plasma ap- 
proach, and is the largest in absolute value of all negative ab- 
sorptions, the OM and XO may be observable. We assumed 
in the previous section that the largest in absolute value 
negative absorption coefficient dominates, and that all the 
energy goes into amplifying this mode. With this assump- 
tion we derived the saturation length and the amplification 
for this dominant mode. This is similar to the assumption 
in previous studies, which give the dominant mode in some 
parameter region as the mode with the highest growth rate 
(White et.al. 1986; Aschwanden 1990a), and it is justified be- 
cause of the exponential nature of the intensity. There are, 
however, two scenarios where modes with smaller absorption 
coeffici ents can be observed. The first scenario, presented in 
section 6.2, depends on the high brightness temperatur e of 
the maser. The second scenario , presented in section 6.3, 



depends on the slowness of the Z-mode. If the ZM is slow 
enough, another mode may be able to reach its full amplifi- 
cation before the ZM saturates, and then, for a short time 
until the ZM saturates, this mode can also be observed. 

Finally, there is a region where the ZM does not have 
negative absorption. For the region 1.5 < v p /vb < 1.6 the 
lowest frequency for which the ZM exists is higher than vb 
and the highest frequency for which the ZM exists is lower 
than 2vb- Negative absorption usually appears near the har- 
monics of the cyclotron frequency, and for the ZM frequen- 
cies with this region of v p /vb there are no negative absorp- 
tion coefficients. 

All the above considerations are summarized in Fig. [j] 
and in Fig. ^ giving the maximum absorption coefficients 
and growth rates for the modes. The coefficients are the 
largest in all frequencies and angles for the given v v /vb, wave 
mode, and with T = 5 • 10 6 K. The ZM does not appear for 
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1.6 < u p /vb because where negative absorption reappears 
at u p /ub > 1.6 the frequencies fail the cold plasma criteria. 
It is important to note that the largest growth rate does not 
necessarily occur at the same angle and frequency as the 
largest in absolute value negative absorption coefficient, as 
can be seen in Tables |l| to |(| 

Whatever mode is amplified, the ECM is difficult to 
observe because of the narrow bandwidth, and especially 
because of the small angular width of the emission. Fig. H 
shows contours of negative absorption coefficients around 
the area where the OM has the largest negative absorption 
(largest in absolute value) . To observe the emission it is nec- 
essary to be within a cone of less than 1 degree around the 
angle of maximum amplification, and within 1 per cent of 
the frequency of maximum amplification. The absorption 1 
per cent away from the maximum drops by half, which for an 
exponential growth with a maximum power of the exponent 
of about 30 means a drop of 7 orders of magnitude in am- 
plification. In Fig. ^ and Fig. |H| the same narrow bandwidth 
and small angular width are seen for the other modes, and 
the phenomena occurs for all other v p /vb values. This gen- 
eral result means that observation of ECM should be rare, 
even if ECM is very common and occurs at all times. 

6.1 Thermal Effect 

In most cases the thermal plasma tends to suppress maser 
emission, as was noted in previous studies (Aschwanden 
1990b; Kuncic & Robinson 1993). For example for v p /vb = 
0.8 (Table |TJ) and T = 5 ■ 10 6 K, the most important OM 
was found at frequency v/vb = 1.03. At a higher temper- 
ature of 10 7 K the most important OM moved to a higher 
frequency of v/ub = 1.079. The fast electrons still gave the 
same negative absorption at v)vb = 1.03, only for T = 10 7 
the thermal absorption was large enough to turn the absorp- 
tion at this frequency positive. On the other hand, the ZM 
at p/vb = 1.2745 was not affected by the thermal absorp- 
tion, because of its distance from the cyclotron frequency. 
The change there was the result of the application of the 
cold plasma criteria, which for T = 10 this frequency vi- 
olates, and therefore the largest acceptable absorption was 
at a lower frequency, further from the resonance and with a 
smaller refraction index. All these changes were not enough 
to change the dominant mode, which remained the ZM, be- 
cause the OM absorption was reduced by a factor of 8, while 
the ZM absorption was reduced only by a factor of 4. The to- 
tal effect of the higher temperature was therefore to reduce 
the relative importance of the OM emission. 

The other ZM given in Table [j] is the ZM wave with 
the largest growth rate. This is a wave at a much lower 
frequency, and a slightly different angle than the wave 
with largest (negative) allowable absorption coefficient. This 
mode also changes with the increase in temperature from 
frequency v/ub = 1.0351 for T = 5 • 10 6 to frequency 
vjvB = 1.087 for T = 10 7 . In this case the absorption 
was only reduced by a factor of 2, and the growth rate by 
about the same. This is an example where the mode with 
the largest growth rate is not the dominant mode. 

A counter example is demonstrated in Table |^. The first 
row of Table |^ gives a negative absorption coefficient for 
the frequency v/ub = 2.04885, which is very close to the 
ZM resonance for this case. The refraction index for this 



frequency and for the angle cos(#) = 0.06 is n re f = 4.42, 
which was large enough for this solution to barely violate 
the 1st plasma criterion (eq. ^), giving a value only 9.5 times 
larger than 0.001683/2. This is not allowed, which is why the 
row is marked as ZMN in the table. Therefore, for v p jvB = 
1.8 the temperature leads us to discard the ZM waves as 
unimportant, and leaves us with the OM as the dominant 
mode. 

When we made the calculations in Table ^ with a tem- 
perature of T = 10 7 K, all the ZM negative coefficients dis- 
appeared entirely, overcome by the thermal absorption, and 
therefore there was no need to even consider the cold plasma 
criteria. 

It is important to note a problem with Table ^ which is 
that the negative absorption coefficients are very small. The 
saturation lengths were therefore very large, and may violate 
the assumption that the maser is very small compared to 
the dimensions of the loop. For example for T — 10 7 K the 
saturation length for the OM was found to be Z sat = 1.8 • 
10 8 cm, which is large enough for our assumption that the 
magnetic field, the density, and the temperature are constant 
to be incorrect. In such a case the maser does not reach its 
saturation intensity, and the loss cone may relax by some 
other means. 



6.2 Non-Dominant Mode 

Since the intensity 7 is difficult to measure we would like to 
estimate the brightness temperature 

(37) 



T b = TTT r~ = 3 - 26 ■ 10 18 . * „ 
2v z Kb (v/GHz) 2 

We can express the intensity in terms of the saturation 
length, use equation|32| and write the intensity (equation |30[) 

as 



j (0.1E k v s \k\ 



AD.A 



v J 



(38) 



Using this, the brightness temperature for the dominant 
mode can now be written as 



Tb 



3.26 • 10 1S 0.1E k v g 



(39) 



(v/GUz) 2 AQAiy 

Therefore for the dominant mode the brightness tempera- 
ture depends only on the mode (through v s ) and on the 
energy assumed to be available to the maser, and not on the 
saturation length Z sa t- 

We estimated the importance of any wave with negative 
absorption by calculating a saturation length assuming this 
particular wave is the only amplified one. We then compared 
the various non-dominant saturation lengths to the satura- 
tion length of the dominant mode, where by definition the 
dominant mode is the mode with the shortest saturation 
length. The amplification length of a non-dominant mode is 
limited to the loss-cone length determined by the dominant 
mode /dom- The intensity of the non-dominant mode can 
then be calculated by using the non-dominant absorption 
and emission coefficients with the dominant mode's satura- 
tion length. The intensity to be used in equation |37| is for 
this case 



T _ Jnon Jfc„ 
1 non — j. I c 



(40) 



© 0000 RAS, MNRAS 000, 000-000 



10 A. Stupp 



or, alternatively 



^non — 



\kn 



0.1E k v g \k 



fe„on/fedo 



(41) 



Where fcdom,Jdom are for the dominant mode, which de- 
termines the saturation length, and fc non ,j non are for some 
other mode, at a different frequency and another emission 
angle, but which also has a negative absorption, and is there- 
fore amplified. This formulation is true no matter which is 
the dominant mode, also for the ZM. 

If |fcnon| ■ idom is large, the resulting intensity is large 
and observable. For example for u p /ub = 1 (Table |^) and 
T = 5 ■ 10 6 K the XO saturation length was the shortest of 
all modes, though very close to the ZM saturation length. 
However, the OM absorption coefficient was also very close 
to the XO absorption coefficient, and obviously the non- 
dominant situation is relevant. This is also one of the rare 
cases where the OM was slower than the ZM, so the growth 
rate of the OM was much smaller than the growth rate of 
the ZM, and certainly smaller than the growth rate of the 
XO. In spite of that, our model predicts a strong emission 
in the OM for this table's u p /vb = 1. For this case j/k for 
the OM is about 1/16 of the j/k ratio for the XO. Therefore 
even when the power of the exponent is the same, the total 
intensity in the OM is only 1/16 of the intensity in the XO. 

If the non-dominant mode has a much smaller absorp- 
tion, for example as in v^jv-g, = 1 (Table ^) and T = 10 7 
K, where /c non ss /cdom/3, then the non-dominant intensity is 
much smaller. For example if we assume fcdom'^dom ~ 25 30, 
as it often is, and the same j/k for the dominant and 
non-dominant modes, the factor between the intensities is 
e 17 ~ 10 7 . This factor reduces a brightness temperature of 
10 18 K to a brightness temperature of only 10 11 K, which is 
lost against the background of the entire loop, and is in any 
case much lower than is observed for microwave spikes. 

The way to decide when a non-dominant mode could be 
observed is to compare absorption coefficients. In Fig. ^|the 
XO growth rate for p^/v-r = 1 (of the example above) ap- 
pears to be much larger than the OM growth rate. However, 
in Fig. |l] we see that the absorption coefficients are equal. 
The example of the previous paragraph shows that the OM 
is observable whenever such a situation occurs. 



6.3 Slow Mode Saturation 

Another possibility for observation of a non-dominant mode 
is for the case of the dominant mode being slow. The Z- 
mode is usually a slow wave, especially near its resonance 
frequency, where its group velocity can be even as low as 
c/500. It is therefore possible that when a certain set of 
conditions is met another mode is amplified to its maximum. 
This can happen if the Z-mode is slow enough that the other 
mode reachs its saturation length before the ZM saturates. 
In that case until the ZM saturates and reduces the enegry 
in the distribution, the other mode is emitted with the full 
amplification of its own saturation length. 

The growth rate equation [l^ already includes the group 
velocity, and therefore a comparison of growth rates suffices 
to show if this slow-mode scenario is possible. It is also pos- 
sible to take the group velocity into account from the per- 
spective of the spatial absorption coefficient, though we need 



a few assumptions. We assumed a typical emitting electron 
energy of 7 = 1.1, and that the important electrons are on 
the edge of the loss cone with pitch angle cos(<^)) = 0.8. These 
parameters gave a velocity parallel to the magnetic field of 
c/3 for the most important electrons, while the OM or XO 
wave was usually faster than c/3. 

For an example we used Table ^| for temperature of 10 7 
K. The OM group velocity is 0.783 • c, and the ZM group 
velocity is 0.1253- c. An OM wave of speed 0.783 c beginning 
near the mirror point at time t\ after the first electrons mir- 
rored overtakes these first electrons (of parallel velocity c/3) 
at time ti = 1.74 ■ t\, and at a po int 0.58 c • t\ from the mir- 
ror point. In our model (section 4.2), when this distance is 
the saturation length, 0.58 • c- 1\ — l°J? , the intensity of the 
radiation is so great that the end-most volume element with 
a loss-cone distribution is drained of its free energy, and the 
loss-cone distribution disappears from that end-most vol- 
ume. A Z-mode wave emitted at time t = 0, and having a 
group velocity of 0.1253 c, has propagated by this time only 
to a distance of 0.1253 ■ c ■ t 2 = 0.376 ■ l°Jf . If the ZM with 
this group velocity has an absorption coefficient such that 
its saturation length is for example 0.4 of the OM satura- 
tion length, the Z-mode had not saturated by the time the 
OM emission has reached its peak and is beginning to close 
the loss-cone. However, at some later time the ZM will have 
had time to saturate near the mirror-point, and has closed 
the loss-cone close to the mirror. An OM/XO wave starting 
from the mirror point at that time does not see a loss-cone, 
and is not amplified. For the case in Table |i| the ZM has 
a saturation length only 0.27 of the OM saturation length, 
and therefore the ZM has saturated by the time the OM 
wave reaches its saturation length. This however does not 
mean that the OM is not amplified as explained below. 

A better calculation than the above takes into account 
that the OM wave advances faster than the electrons and 
much faster than the ZM wave. The requirement then is that 
the electrons the OM wave passes have not themselves gone 
previously through a saturated Z-mode wave, and therefore 
have retained the loss-cone distribution. The result of that 
requirement is that for a mode with i sat to be amplified, the 
Z-mode saturation length Z z must be large enough so that 
the inequality below is correct : 



h > V z V g — Ve 
Zsat ~ V e U g — V z 



(42) 



With v e the velocity of the electrons, v s the group velocity 
of the OM or XO, v z the group velocity of the ZM, Z z the 
saturation length of the ZM, Z sa t the saturation length of the 
OM. The same equation holds for any two modes, one with 
group velocity v z slower than the electrons, and the other 
with group velocity v s faster than the electrons. 

For the case of Table ^ the left side of the equation 
is 0.27 while the right side is 0.257, and therefore the OM 
can reach its maximum possible amplification before the ZM 
closes the loss-cone and prevents it. 

A simpler estimate for the relative importance of the 
modes is to compute r = Z sat /« g for each of the modes. 
This estimate is very close to calculating the growth rate, 
since dividing the saturation length by v s is equivalent to 
replacing the absorption k by the growth rate F = —k - v g in 
equation |x| This criterion tends to be smaller for a larger 
growth rate, and the smallest value is the most important. 
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For our example (Table ^, T = 10 7 ) the 7th column of the 
table shows that the growth rate of the OM is larger than 
the growth rate of the ZM. Therefore, even though the sat- 
uration length of the OM is much larger than the saturation 
length of the ZM, and after the ZM saturates the OM disap- 
pears, the growth rate criterion, and our calculations above, 
show that the OM can be observed for a short duration, 
until the ZM saturates . 

Another example is for v p /vb = 0.8 (Table jjj), T = 
5 • 10 6 and the ZM. In the Table we give both the ZM wave 
with largest (absolute value) negative absorption, and the 
ZM wave with the largest growth rate. The ZM wave with 
the largest absorption is very slow and has a very small 
growth rate, and therefore the important wave is at first the 
one with the largest growth-rate. For v p /vb = 0.8 the largest 
growth rate in the ZM is still smaller than the growth rate 
of the OM wave, and therefore we conclude that the OM 
reaches its saturation length before the ZM closes the loss- 
cone, and we expect to see an OM wave for a short time. It 
turns out that the difference between the time to saturation 
for the OM and the time to saturation for the largest growth 
rate ZM is very small, about 0.1 msec. The time we expect 
to see the OM at its full amplification is therefore very small. 

In this case it was usefull to check the non-dominant- 
mode approach of the previous section, and we found that 
the OM is observed until the small saturation length ZM 
saturates. When we used the OM absorption coefficient with 
the saturation length of the largest growth rate ZM we got 
a brightness temperatures of order Tb = 10 15 , which is ob- 
servable. 



6.4 Time Dependence 

The duration of the pulse of fast electrons is obviously un- 
related to the loss-cone time scales. If the pulse is shorter 
than the loss-cone time-scale, less than t 2 = 3 • Z sa t/c for 
the example of the previous section, the maser intensity in 
the OM is not strong enough to close the loss-cone. The 
loss-cone distribution therefore moves from the mirror point 
at the typical loss-cone velocity, and the maser continues 
to emit, though at a weak intensity. Using the terminology 
of the example from the previous section, and assuming a 
pulse duration of ti, the maximum amplification length is 
c-ti/3. In the example we had Z sa t = 0.58-c-ti , and therefore 
the maximum amplification length is about Z sa t/2. Since the 
intensity is exponential with the amplification length, this 
translates into a much weaker emission. Typically the power 
of the exponent for the dominant mode is 20-30, and with 
Zsat/2 the intensity is more than a hundred thousand times 
weaker. For such a case it is unlikely that the maser emission 
can be observed at all. 

The other extreme is that the pulse of fast electrons per- 
sists for a time much longer than typical maser saturation 
times. For such a case, the loss-cone closes at distance Z aa t 
from the mirror point, but new fast electrons move into place 
constantly. At first glance it would seem that the maser per- 
sists for approximately the duration of the pulse, however, 
if there is a ZM with a smaller saturation length than the 
observable modes this does not happen. After the Z-mode 
reaches its saturation length it quenches the fast electron 
loss-cone distribution, and the observable modes can not be 
amplified enough to be seen. In this picture the ZM waves 



saturate and maintain an intensity just strong enough to 
drain every new volume element of its energy when it over- 
takes the ZM waves. An observable OM or XO wave starting 
from the mirror at a time t > £i of the previous section, so 
that the ZM has already reached saturation t > Z^t / v g M > 
does not see a loss-cone, and therefore is not amplified. For 
our example above, the intensity of an observable mode at 
first increases, reaching a maximum at time £2 after the first 
electrons are mirrored. At this time the OM wave which left 
the mirror at t\ reaches Z sat . This high intensity can last 
only until OM waves leaving the mirror encounter saturated 
ZM waves. For the example from Table ^ the ZM reaches 
its saturation length near the mirror at time 1.25 • t\. An 
OM wave leaving the mirror point at 1.25 • t\ reaches Z sa t 
at t% + 0.25ti, and the maser then starts to decrease in in- 
tensity. The maser is therefore at peak intensity only for a 
time 0.25ii = 1.7 msec. Of course, if the shortest saturation 
length is for an observable mode, the maser does persist for 
the length of the pulse of accelerated particles. 

Our conclusion is therefore that if the smallest Z sat is 
for the ZM, the maser persists for a time less then but of 
order [l aa ,t/v s ]zM, regardless of the duration of the fast elec- 
tron acceleration pulse. This is the case for Table [j] at all 
temperatures. Of course, after the pulse ends, the waves dis- 
appear, and when a new pulse reaches the mirror point, the 
process can start again. 

An example of a different kind is for v^jv-e, = 1 (Table 
^ and T = 5 ■ 10 6 K. For this case the saturation length of 
the XO was shorter than the saturation length of the ZM, 
and the saturation length of the OM was not much longer. 
Therefore we expect the XO radiation to last, and the OM 
radiation to be amplified as well. For the other tables, at 
T = 5 • 10 6 the smallest saturation length was always for an 
observable mode, and therefore at this temperature the loss 
cone should always produce long lasting maser emi ssion . 

The differential approach described in section |4.2| also 
gives a model for the time-profile of the maser. The emis- 
sion should appear as a rising exponent, and then decrease 
faster than an exponent, because the elements whose energy 
was exhausted begin to positively absorb radiation from the 
still amplifying elements. This picture is at odds with the 
characteristics observed, which were a rise a s e'i 2 ) and a 
decrease as an exponent (Sudel & Benz 1990). Perhaps the 
gaussian rise is a feature of the acceleration mechanism, for 
example if rihot also increases with time. 



7 COMPARISON WITH PREVIOUS WORK 



We can compare our results with only one previous case ( As 
chwanden 1990b) where saturation lengths were estimated 
and given for example for v v jua — 0.5 as ifJi' — 10 km and 
Z s ° a f = 30 km. 

We incorporated in our program Aschwanden's thermal 
distribution in energy and sin 6 (3 x 0) loss-cone. Aschwanden 
used the standard parameters Tkot = 10 8 K, Tcoid = 10 6 K, 
JT-coid = 1-25 • 10 8 cm -3 , and n hot /n co id = 0.01. With the 
above and from the standard v-p/us = 0.1 a magnetic field 
of B = 360 gauss was the standard field. If the number of 
particles was kept constant the magnetic field for v^jv-e, = 
0.5 as in the example above should have been reduced to 
B = 72 gauss. 
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Our calculations with Aschwanden's distribution func- 
tions gave the normalized results, with B — 100 gauss 
and divided by rihot, for the OM: k = —1.1 ■ 10 -11 , j = 
1.69 ■ 10" 24 , c/vg = 1.144, u/u B = 1.019, cos(0) = 0.24 and 
for the ZM: k = -3.6 • 10" 11 , j = 4 ■ 10" 23 , c/v s = 7.3, 
u/i/b ~ 1.039, cos(#) = 0.2. Where these numbers are for 
the respective largest growth rate in every mode, since that 
was Aschwanden's criterion for the dominant wave. 

Using our simple method to calculate the saturation 
length, and the scaling laws k oc nh ot /B and j oc rihot x B, 
we got l^f = 4.89 km and 1°^' = 17.9 km. These saturation 
lengths are about a factor 2 smaller than Aschwanden's. We 
also have to take into account Aschwanden's assumption of 
an initial energy density equivalent to Tb = 10 14 K. In our 
model the waves amplify the spontaneous emission, and the 
amplification is according to the radiative transfer equation 
|3c| . A typical value is j/k « 10 -13 and if this is the intensity 
I the equation for the brightness temperature (equation |37[ ) 
gives the equivalent temperature as only Tb = 10 7 K. The 
brightness temperature reaches Tb = 10 14 K only after the 
wave has been amplified by a factor of 10 7 w e 16 , and the 
comparison of the saturation lengths should be only from 
this stage. 

The total saturation length gives a power in the expo- 
nent of about 30, therefore if we had started with a wave 
at Tb = 10 14 K the saturation length would have been only 
half of what we calculate. The conclusion is that computing 
the saturation length with our method resulted in satura- 
tion lengths of a factor of 4 smaller than those calculated 
by Aschwanden's method. This was the result of neglecting 
the change in the loss-cone during the amplification of the 
wave, and assuming that the absorption coefficient remains 
constant until the wave reaches its maximum energy. The 
differential iterative method shows that the actual satura- 
tion length is up to twice what the simple method gives, 
and therefore we get a factor of about 2 back. Considering 
the simplicity of the equations we use, a factor of two is 
extremely close. 



8 SUMMARY 

We have developed a model for estimating the intensity and 
duration of emission through the Electron-Cyclotron-Maser 
(ECM) mechanism. 

Observations led us to use a power law for the fast elec- 
trons, with a minimum energy of 10 keV and maximum en- 
ergy of 1 MeV. This distribution is harder than the ones 
studied in the past (Winglee & Dulk 1986; Aschwanden 
1990a). 



From recent results (Aschwanden & Benz 1997) we esti- 
mate the important ratio v p /pb in the foot-points of flares, 
which is the region of interest for this study, to be in the 
range 0.5-3. For most of this range a loss cone in a power law 
distribution of electrons results in the ZM being the domi- 
nant magneto-ionic mode. The ZM is unobservable, because 
when a ZM wave reaches a region with a smaller plasma 
frequency than the plasma frequency where the wave origi- 
nated, it is absorbed. Thus it seems at first glance that the 
maser mechanism can not produce observable emission. 

We use an absorption coefficient approach, differing 
from previous studies which used the growth rate. Using this 



approach we calculate the saturation length for the compet- 
ing magneto-ionic modes, and show that in many cases the 
observable OM and XO are amplified to observable intensi- 
ties. 

The conditions which enable the OM and XO modes 
to be amplified in the presence of a ZM with much larger 
(in absolute magnitude) negative absorption are : first, the 
inclusion of an ambient thermal plasma with temperature 
of order T = 10 7 K ; second, the low group velocity of the 
ZM which enables a competing mode to reach its saturation 
before the ZM saturates ; third, the high amplification of the 
maser, which makes even non-dominant modes observable. 

Our conclusion is that in the range v p /vb = 0.5 — 1.8 
the maser mechanism can produce observable emission in 
the form of intense spikes with a narrow bandwidth - about 
one per cent of the frequency - and narrow angular width - 
about 0.17T solid angle. 

The examples given in the tables indicate that long du- 
ration emission in the OM and XO appears at angles of 30- 
60 degrees to the magnetic field, and at 3-6 per cent of the 
cyclotron frequency above the nearest cyclotron harmonic. 
This is usually the second harmonic for the OM, and always 
the second harmonic for the XO. We expect the long dura- 
tion emission to be determined by the duration of the pulse 
of accelerated electrons, and the brightness temperature to 
be of order Tb = 10 16 K and above. 

Short duration emission appears when the OM or XO 
are emitted only until the ZM saturates. For these cases the 
saturation lengths of the OM and XO are longer than the 
ZM saturation length, but the ZM group velocity is smaller 
than the group velocity of the OM or XO. The short duration 
emission lasts in our model from fractions of a millisecond 
to a millisecond, for example for v p /i/b = 0.8 (Table [l]) and 
a temperature of 5 • 10 6 K, where the OM has the largest 
growth rate. 

From the results given in Tables | to | for a temperature 
of 5 ■ 10 6 K, we expect emission for the parameters of the 
tables for the duration of the accelerated electrons' pulse, 
and at frequency v w 2.06^b- The exception is the case 
of Pp/vB = 1 (Table |^) where the XO emission for T = 
5 ■ 10 6 K is at u/vb ~ 2.06, and there is also OM emission 
of observable intensity at v/vb ~ 1.03. We also expect long 
duration emission for v p /ub = 1.2 (Table and T — 10 7 



K. For this case U 



and therefore the XO will 



be amplified to observable intensities. The XO frequency is 
v = 2.129^b- For the cases v p /vb = 1.6 and v p /vb = 1.8 
(Table | and Table |) and T = 10 7 and T = 2 ■ 10 7 , the OM 
is the dominant mode, but the saturation lengths are much 
too large. We conclude that the loss-cone relaxes by some 
other means, and not through the ECM mechanism. 

In summary, this study shows that the ECM produces 
millisecond radio spikes such as observed from the Sun. The 
ECM emission reproduces the frequency of the observed 
spikes, the band width of the spi kes, the source size as re- 
cently observed (Benz et.al. 1996), the brightness tempera- 
ture, and the duration of the spikes. 

The ECM does not produce emission at two frequen- 
cies simultaneou sly, and therefore can not reproduce the 
observed bands (Krucker fc Benz 1994). However, if a loss 



cone distribution forms independently in different regions, 
and the magnetic field is not the same in these regions, then 
ECM in several frequencies and angles can be produced. For 
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example in Table |2J for temperature T = 5- 10 6 the XO emis- 
sion is at angle 69 degrees. In Table ^| for T = 5 ■ 10 6 the 
XO emission is at angle 64 degrees. If two regions, one with 
Vp/vs = 1 and one with u v /vb = 1.2, emit ECM, and if the 
magnetic field changes direction by a few degrees between 
the regions, it may become possible to observe both. Since 
the magnetic field and the density both increase towards the 
foot points, the ratio between the emitted frequencies in the 
two regions can be larger than 1.2. 

Our mode l also does not repr oduce the reported gaus- 



sian rise time ( Gtudel fc Benz 1990 ). However, if the number 



of fast electrons also increases with time, the rise time is a 
feature of the acceleration mechanism, since the emission in 
ECM is directly proportional to the number density of the 
fast electrons. 

It is possible to conclude from our results that the ob- 
served frequency of the spikes should be higher at the begin- 
ning of a flare, since the hot and dense chromospheric plasma 
is then still close to the chromosphere, where a large mag- 
netic field causes a large cyclotron frequency. As the flare 
advances, and the evaporation front propagates upwards, 
temperatures below increase, the unobservable ZM domi- 
nates for the larger temperatures, and the magnetic field in 
the region from which ECM emission is observed is smaller. 



© 0000 RAS, MNRAS 000, 000-000 



14 A. Stupp 



REFERENCES 

Aschwanden M.J, 1990a, A&AS, 85, 1141 
Aschwanden M.J., 1990b, A&A, 237,512 
Aschwanden M.J., Bcnz A.O. , 1988, A&A, 332,447 
Aschwanden M.J., Bcnz A.O. , 1997, ApJ, 480,825 
Bcnz A.O., 1986, Sol. Phys., 104, 99 

Benz A.O., 1993, Plasma Astrophysics Kinetic Processes in Solar 

and Stellar Coronae. Kluwer Academic Publishers 
Benz A.O., Graham D., Isliker H., Andersson C, Koehnlein W., 

Mantovani F., Umana G., 1996, A&A 305, 970 
Fleishman G.D., Yastrcbov S.G., 1994a, Sol. Phys., 154, 361 
Fleishman G.D., Yastrcbov S.G., 1994b, Sol. Phys., 153, 389 
Gudcl M., Benz A.O., 1990.A&A, 231,202 
Holman G.D., Benka S.G., 1992,ApJ, 391, 854 
Holman G.D., Eichler D., Kundu M.R., 1980, in M.R. Kundu, 

T.E. Gergely , eds, Proc. IAU Symp. 86, Radio Physics of the 

Sun. Dordrecht : Reidel, p. 457 
Kuncic Z., Robinson P.A., 1993, Sol. Phys., 145,317 
Krucker S., Benz A.O., 1994,A&A, 285,1038 

Lang K.R., 1980, Astrophysical Formulae 2nd Edition. Springer 
Verlag 

Mackinnon A. , Vlahos L., Vilmcr N., 1992.A&A, 256, 2, 613 
Melrose D.B, 1989, Instabilities in Space and Laboratory Plasma. 

Cambridge University Press 
Melrose D.B, Dulk G. A., 1982, ApJ, 259, 844 
Ramaty R., 1969, ApJ , 158, 753 
Trulsen J., Fejer J.A., 1970, J. Plas. Phys., 4,4, 825 
White S.M, Melrose D.B., Dulk G.A., 1986, Adv. Space Res. 

6,6,163 

Winglee R.M.,1985, ApJ, 291,160 

Winglee R.M, Dulk G.A., 1986,Sol. Phys., 104,93 

Wu C.S., Lee L.C, 1979, ApJ, 230,621 

Zirin H., 1989, Astrophysics of the Sun. Cambridge University 
Press 



© 0000 RAS, MNRAS 000, 000-000 



Electron- Cyclotron-M 'as er Observable Modes 



Table 1. v p /v B = 0.8 
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Table 6. u p /u B = 1.8 
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r 


Mode 


Kelvin 




Degrees 


cm 2 


cm 


cms -1 


cm 3 s _1 




5 


2.04885 


86.6 


-2.65E-08 


0.00019 


6.45E+07 


1.714 


ZMN 


5 


2.0596 


57 


-1.5E-13 


48.9 


1.44E+10 


2.2E-03 


OM 


5 


3.0844 


68 


-5.4E-14 


122.3 


2.3E+10 


1.1E-03 


XO 


10 


2.1237 


48 


-3.97E-14 


180.8 


1.68E+10 


6.7E-04 


OM 


10 


3.1539 


60 


-1.94E-14 


331.3 


2.08E+10 


4.03E-04 


XO 


20 


2.2893 


34 


-2.24E-15 


3045 


2.18E+10 


4.9E-05 


OM 
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Figure 1. Largest in absolute value negative absorption coeffi- 
cient as function of v p /vb, for the ZM, OM and XO mode. The 
calculation was performed with temperature of 5 • 10 6 Kelvin, 
™cold/ n hot = 10 4 , power-law index of S = 3, low-energy cutoff 
of 10 kcV, high energy cutoff of 1 MeV, loss cone parameters 
cos(«) = 0.81,cos(a — Aa) = 0.83. The absorption coefficient 
given for any v p /vb is the largest in the entire frequency range 
v = z/b — 3vb, and angle to magnetic field range 8 = — 90 
degrees. 




Figure 2. Largest growth rate as function of v p /isb, for the ZM, 
OM and XO mode. The calculation was performed with tempera- 
ture of 5-10 6 Kelvin, n co id/ n hot = 10 4 , power-law index of S = 3, 
low-energy cutoff of 10 keV, high energy cutoff of 1 MeV, loss cone 
parameters cos(«) = 0.81, cos(a — Aa) = 0.83. The growth rate 
given for any fp/fB is the largest in the entire frequency range 
v = v-q — Sz-'b , and angle to magnetic field range 9 = — 90 de- 
grees. The largest growth rate does not necessarily occur at the 
same angle and frequency as the largest absolute value negative 
absorption. 



Figure 3. Contour map of the area surrounding the maximum 
negative absorption in the OM for v p /vb = 1, temperature T = 
5 ■ 10 6 Kelvin, and n cold /n hot = 10 4 . The x-axis of the figure 
denote cos(S), the y-axis denote the normalized frequency v/v-g., 
and the absorption coefficients are normalized to n hot = 1, and 
multiplied by -10 12 
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Figure 4. Contour map of the area surrounding the maximum 
negative absorption in the XO mode for v p /v-g. = 1, temperature 
T = 5 ■ 10 6 Kelvin, and n co id/"hot = 10 4 - The x-axis of the figure 
denote cos(0), the y-axis denote the normalized frequency v/v-g., 
and the absorption coefficients are normalized to n hot = 1, and 
multiplied by -10 12 
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Figure 5. Contour map of the area surrounding the maximum 
negative absorption in the ZM for v p /v-g, = 1, temperature T = 
5 ■ 10 6 Kelvin, and ri cold /ra hot = 10 4 . The x-axis of the figure 
denote cos(#), the y-axis denote the normalized frequency u/u-g, 
and the absorption coefficients are normalized to n hot = 1, and 
multiplied by -10 12 
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